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Abstract 

We present results of a detailed study of the three-body T,c{24:55)^pTT~^ decay. 

A significant enhancement of signal events is observed in the p7r~^ mass system 
near 1.5GeV/c^ that is consistent with the presence of an intermediate baryonic 
resonance N^, where iV° is the 7V(1440)°Pii or iV'(1535)''S'ii state, or an admixture 
of the two states. We measure the product Sc(2455)°iV°) x B{N^ 

p7r+)=(0.80 ± 0.15(stat.) ± 0.14(syst.) ± 0.21) x 10"^, where the last error is due to 
the uncertainty in B(A'^ pK~TT~^). The significance of the signal is 6.1 standard 
deviations. This analysis is based on a data sample of 357 fb~^, accumulated at 
the T(45) resonance with the Belle detector at the KEKB asymmetric-energy e'^e~ 
collider. 

Key words: B-meson, Charmed baryon, Pu and resonaces. 

PACS: 13.20.H 



* Corresponding author. 

Email addresses: hokimSknu.ac.kr (H. O. Kim), kichinii@post.kek.jp 
(H. Kichimi). 



4 



Various charmed baryonic B decays into four-, three- and two-body final states 
have been reported p!|2|3Pf5|6f7] , and the measured branching fractions show 
clearly that the branching fraction increases with the multiplicity of the final 
state [HfS]. To understand this hierarchy, it is interesting to study decays of 
pvT^TT into three- and two-body final states. The branching fractions 
are predicted from CKM matrix elements [9], while the form factors of the 
decay vertices depend on the decay mechanism. Experimental studies provide 
stringent constraints on the theoretical models p^ll|12] . 

In this report, we perform a detailed study of the intermediate three-body de- 
cay B^ — s> Sc(2455)''p7r^ observed in the previous analysis of B^ A^pvr'^vr" pQ, 
using a data sample of 388 x 10^ BB events, corresponding to 357 fb^^ accumu- 
lated at the T(45') resonance with the Belle detector at the KEKB asymmetric- 
energy e~^e~ collider [T3] . 

The Belle detector is a large-solid- angle spectrometer based on a 1.5 Tesla su- 
perconducting solenoid magnet. It consists of a silicon vertex detector (SVD) 
(a three-layer SVD for the first sample of (152.0 ± 1.2) x 10^ BB events 
and a four-layer SVD for the latter (235.8 ± 3.6) x 10^ BB events), a 50- 
layer central drift chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel-like arrangement of time-of-fiight scintillation coun- 
ters (TOF), and an electromagnetic calorimeter (ECL) comprised of Csl (Tl) 
crystals located inside the superconducting solenoid coil. An iron fiux return 
located outside the coil is instrumented to detect mesons and to identify 
muons (KLM). The detector is described in detail elsewhere [H]. We simulate 
the detector response and estimate the efficiency for signal reconstruction by 
Monte Carlo simulation (MC). We use the EvtGen program [1^ for signal 
event generation and a GEANT-based [16] detector simulation program to 
model the Belle detector response for the signal. 

We first describe briefiy the previous analysis of B^ — > A+pvr+vr" [T], We select 
B^ —>■ A'^p7i~^7r~ events by reconstructing A+ —>■ pK~TT^ decays, using charged 
tracks reconstructed by the SVD and CDC, and hadron identification informa- 
tion (such as protons, kaons and pions) provided from the CDC dE/dx, TOF 
and ACC (PID) [17], and ECL and KLM information to veto electron and /i 
tracks. Charge-conjugate modes are implicitly included throughout this paper 
unless noted otherwise. After the event selection, we fit the AE distribution for 
the B candidate events with 5.27 GeV/c^ < Mbc < 5.29 GeV/c^, with a double 
Gaussian fixed to the signal MC shape (acoro = 7MeV/c^, o"taii = 16MeV/c^) 
plus a linear background. The variable AE = Eb — -Ebeam is the difference be- 
tween the reconstructed B meson energy (Eb) and the beam energy (-Ebeam) 
evaluated in the center-of-mass system (CMS), while Mbc = \/ -^bcam ~ 
is the beam-energy-constrained B meson mass and Pb is the momentum of 
the B meson also evaluated in the CMS. We obtain a B signal of 1400 ± 49 
events for B^ — > A+p7r+7r~. Figure [1] shows the A+tt"'' and A+tt" mass dis- 
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Fig. 1. The mass distributions of (a) A^tt"*" and (b) A^vr" in —>■ A^pvr^Tr". 
The points with error bars show the mass distribution for the events in the B signal 
box, and the shaded histogram indicates that for the background. The sohd and 
dashed curves represent the signal and the background, respectively, obtained from 
a simultaneous binned likelihood fit. 



tributions for the events in the B signal region \ AE\ < 0.03 GeV (±4(t) and 
5.27GeV/c^ < M^c < 5.29GeV/c^ (±4o-). We focus our discussion on the 
Ec(2455)+^ and Sc(2455)° resonances clearly observed in figure [Ij; (182ibl5) 
events for 5° ^ Ec(2455)++p7r- and (122±14) events for 5° ^ Sc(2455)°p7r+, 
corresponding to branching fractions of (2.1±0.2(stat.)±0.3(syst.)±0.5) xlO~^ 
and (1.4±0.2(stat.)±0.2(syst.)±0.4)xl0~^, respectively. Hereafter, we denote 
Se(2455) as S^. 

Figure [2] shows (a) the Dalitz plot and (b) the M'^{pTX^) distribution for the 
S++j)7r^ events, and (c) the Dalitz plot and (d) the M^(p7r+) distribu- 
tion for the B^ —>■ S°p7r+ events. Here we require the Sc candidates satisfy the 
invariant mass requirement 2.447 GeV/c^ < M(A+7r=^) < 2.461 GeV/c^ (±2a). 
We find that the M^(p7r~) distribution for B^ —>■ S+'''pvr~ is consistent with 
three-body phase space, while the M'^{piT^) distribution for B^ — >■ S^pvr"^ 
has a significant peak. In what follows, we present a detailed study of the 
B^ — s> S°p7r+ decay. 

Figure [3] shows the AE distribution for the B^ —>■ S^p7r+ events, which are 
selected from the B^ —>■ A+pTr+vr" sample with the additional requirement 
that the A^tt^ mass be consistent with the The curves show fits to the 
data with a double Gaussian function with shape parameters fixed to the 
values from signal MC and a linear background. We obtain a B signal yield of 
(102 ± 11) events and a background of (17 ± 3) events. The signal reduction 
of 16% is consistent with the MC estimation of the effect due to the S° mass 
requirement. We estimate a non-E° background of (8 ± 4) events from a fit 
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to the AE distribution in the Sc mass sideband 2.435 GeV/c^ < M(A+7r") < 
2.442 GeV/c^ and 2.466 GeV/c^ < M(A+7r-) < 2.473 GeV/c^. This can be 
compared with (2.5 ± 0.5) events estimated from MC simulation of — >■ 
A+pTT+vr^ decay with four-body phase space normahzed to the total of 1400 
events [1]. Here the error is due to the statistics of the simulation. We do not 
take into account the non-S|? background in the analysis that follows. 
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Fig. 2. (a) Dalitz plot and (b) M^ipTT') distribution for S++p7r". (c) Dalitz 

plot and (d) M'^{pTT~^) distribution for — > S^pvr"'". Points with error bars indicate 
the data, and histograms are the decays simulated according to three-body phase 
space. 
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Fig. 3. AE distribution for the B^ — > A^pvr^vr^ events in the signal region 
with 2.447 GeV/c^ < M(A+7r=^) < 2.461 GeV/c^. The curves indicate the fit with a 
double Gaussian for the signal and a linear background. 

Figure H] shows (a) the pvr"*" mass, (b) cos 9p and (c) ScTt^ mass distributions 
for the selected S°p7r+ events. Here, cos6'p is the cosine of the angle 

between the p momentum and the direction opposite to the B momentum in 
the pvr"'" rest frame. The shaded histograms indicate the distributions for the 
background discussed above. The background shapes are obtained by fits to 
the data in the sideband region \AE\ < 0.1 GeV and 5.26GeV/c^ < Mbc < 
5.29GeV/c^ outside the B signal region, and the yield is fixed to 17 events. 
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Here, the M{p7i~^) distribution is parameterized by the function Pbkg(^) = 

C ■ A/tmin ■ Uax " (1 + Ci ■ tmin)(l + C2 ■ tmax) with t^in = (M^ - M^^J and 

tma.x = (^max~^^)- ^min and Mmax are the minimum and maximum masses. 
The variable c is a normahzation constant, and Ci and C2 are shape parameters. 
The cos 9p distribution is modeled by a second-order Chebyshev polynomial. 




MOJii*) (GeV/c? ) 



Mp:.*) (GeV/<?) 



Fig. 4. Data distributions for (a) M(p7r+), (b) cos6'p and (c) M(S^7r+). The shaded 
histograms indicate the normaUzed background. 

We find a significant structure in the pvT"*" mass distribution, and a forward 
peak in the cos 6p distribution, and a low mass ScVr"*" enhancement, denoted 
as (S°7r"'")x. The pvT"*" mass structure has a mass near 1.5 GeV/c^ and a width 
of about 0.3 GeV. We denote it as iV°, and investigate its characteristics in 
detail. In order to describe the pvT"*" mass structure, which is not explained 
by a simple phase space non-resonant H^pn^ decay, we consider an 

intermediate two-body decay S° — > S^iV" with a resonant state iV° pix^ . 
However, we still cannot reproduce the forward cos 6p peak and the S°7r+ low 
mass structure with these two modes only. Therefore, we introduce one addi- 
tional mode (S°7r"*")xp to account for the observed features. As the low 
mass Sc7r+ structure is close to threshold, it produces a forward peak in the 
cos^p distribution. In the low E^tt"*" mass region, we search for known S°7r+ 
resonant states ^rE\ in finer mass bins, but find no signals. So far, there is 
no good candidate to interpret this broad structure as a resonance. There- 
fore we assume that there is a threshold mass enhancement with a mass of 
2800MeV/c2 and a width of 350 MeV obtained from a fit to the E°7r+ mass 
distribution using a relativistic Breit-Wigner (S-wave) function. 

Figure compares binned Probability Density Functions (PDF) of the MC 
simulated events for the three assumed decay modes. The histograms show the 
PDFs for (a) the M(p7r+), (b) cos^p and (c) the M(S°7r+) distributions. The 
solid histograms show the distributions of the mode B^ H^N^, assuming a 
P-wave relativistic Breit-Wigner amplitude with a mass of 1530MeV/c^ and 
a width of 340 MeV. 



To determine the iV° mass, width and the yields of the three modes, we per- 
form a maximum likelihood fit to the observed M{p7^~^) and cos 6p distribu- 
tions shown in Fig. |H These two distributions are sufficient to fully describe 
the three-body decay B^ S^pvr"*". To model the observed distribution, we 
construct a function F {M{p'n'~^), cos 6 p) from the sum of PDFs of the three 
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Fig. 5. Binned probability distributions of (a) M{p7r^), (b) cos 9p and (c) M(S[!7r+), 
where we compare MC simulated distributions for — > Ti^N^ (solid lines), 
—>■ S|?p7r+ (dashed lines), and B^ —>■ (S|?7r"'")xp (dotted lines). We make a 
simultaneous fit to the distributions in (a) and (b). 

decay modes and the background. 



where Pi denotes a product of the normahzed PDFs, Qi(M(p7r+)) (20 bins) 
and Ri{cos6p) (16 bins), and z/j stands for the yield of the z-th mode. We 
plot Q(M(p7r+)) and R{cos6p) distributions from the detector MC simula- 
tion for — > (S°7r"*")xP and S°p7r+ modes, as shown in Figs [5t^a) 
and (b), respectively. For the B^ S°iV° mode, we use the MC simulated 
R{cos6p) distribution and a Breit-Wigner for Q{M{p7i^)) with the mass 
and width {m^, F) as free parameters. A small systematic error due to the 
use of the Q{M{p'n'~^)) distribution without the MC detector simulation is 
discussed later. We use a P-wave (S-wave) relativistic Breit-Wigner shape. 
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Here m is the mass of the pvT"*" system, and is the nominal iV" mass, 
and F(m) is the width. The variable p is the momentum of a daughter par- 
ticle in the iV" rest frame, and po is that for the nominal iV° mass. B{p) = 
1/ \Jl + {RpY is the Blatt-Weisskopf form factor [19]. The value i?, called the 
centrifugal barrier penetration factor, is set to 3 (GeV/c)~^ for a P-wave, and 
is zero for an S-wave. L indicates the orbital angular momentum. For the S- 
wave Breit-Wigner amplitude pO] we use F(m) with m = m^j and L = to 
parameterize the smooth shape near the mass threshold. 



We define an extended unbinned likelihood with coarse bins, and carry out a 
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maximum likelihood fit. 



g-{i^l+i^2+!^3+!^4) 



J|F(mR, r, ui, 1/2, z/3, 1/4) 



We fit the mass and width (mij,r) and the yields z/i,i/2 and 1/3 as free 
parameters, while the background yield is fixed to 17 events. Table [1] sum- 
marizes the fit results with various model assumptions. We calculate the statis- 
tical significance from the quality -21n(Lo/-Lmax), where Lmax is the maximum 
likelihood returned from the fit, and Lq is the likelihood with the signal yield 
fixed to zero, and taking into account the reduction of the degrees of freedom. 
We obtain a significance of 7.0 cr for the —>■ S^iV" contribution. The signal 
in the mode (S°7r+)xp has a statistical significance of 4.6 cr, while that 

for B^ — >■ S^pvT"*" is not significant (0.8 cr). Here, we calculate the goodness-of- 
fit from the likelihood ratio A HI 



36 

~ -21nA = 2 ^[J^j(m^, T, u^, z/g, z/3, z/4) - Fj 

+FM^, — ^ v)]' 

JF,(m^,r,z/i,z/2,z/3,z/4) 

where Fj and J^j are the observed and the fitted yields, respectively, in the 
j-th bin: j = 1,20 for 20 bins in Q{M{piT~^)) and j = 21,36 for 16 bins in 
R{cos9p). 

The small contribution (z/3 = — 11 ± 10) can be understood from Fig. [51 
The mode B^ —>■ S^pvT^ has a broad prr^ mass distribution similar to B^ 
(E°7r"'")xP, while it does not reproduce the forward cos^p peak. On the other 
hand, the B° S°iV° mode can reproduce the pn~^ mass bump structure and 
the uniform cos 6p distribution. Hence, we fix z/3 = in the subsequent fit and 
the uncertainty of this contribution is taken into account in the systematic 
error. 
Table 1 

Summary of the simultaneous fits to the M{pTT~^) and cos^p distributions with 
the three decay modes S° Y.°N°, B° ScPvr+ and B° {Y.°Tr+)xp. (a) - 
(d) represent fits with various assumed contributions. Here, we show the fit results 
with the P-wave assumption, as we find no significant difference from the S-wave 
assumption. 



Decay mode 


(a) 


(b) 


(c) 


(d) 


Signif. 




free 


free 


free 





7.0 


50 ^ (S^+)xP 


free 


free 





free 


4.6 


B^ S0p7r+ 


free 





free 


free 


0.8 


xVndf 


31.7/31 


32.4/32 


52.8/32 


88.4/34 





Figure [6] shows the results of a fit to (a) the M{p7i~^) and (b) cos 9p distributions 
under the P-wave assumption. The data are the points with error bars. The 
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Fig. 6. Simultaneous fit to (a) M{pTT~^) and (b) cosOp distributions with a P-wave 
Breit-Wigner. The points with error bars are the data, and the curves are the con- 
tributions from Ti^N^ (dashed), — > (S[!7r"'")xP (dotted), the background 
(shaded) and their sum (solid), (c) M{T,^tt^) distribution, where the curves repre- 
sent their contributions obtained by the fit to (a) and (b). 

curves are the contributions from —>■ S^iV" (dashed), —>■ (S[?7r^)xp 
(dotted), the background (shaded) and their sum (solid). We obtain yields 
of (70 ± 11) and (32 ± 9) for the modes 5° ^ SOiVO and 5° ^ {E'^^ti+)xP, 
respectively. Figure [6t^c) shows that the M(S°7r"'") distribution is consistently 
represented by the fitted parameters even though the M(S°7r"'") distribution 
is not included in the fit. 
Table 2 

The fitted mass and width with relativistic S-wave and P-wave Breit-Wigners. 
The first errors are statistical and the second are systematic including the uncer- 
tainties in the yields of B^ T,^p7r~^ and the background, and the assumption of a 
low mass (S[!7r+)x structure. 



Item 


Yield 


Mass 


r 


xVndf 




Events 


MeV/c2 


MeV 




S-wave 


71 ± 11 ± 10 


1473 ± 31 ± 2 


315 ± 72 ± 53 


32.9/32 


P-wave 


70± 11 ± 10 


1516 ±29 ±14 


365 ± 97 ± 90 


32.4/32 


S^Jvr^ sys. 


±5 


±8 


±50 





Table [2] compares the fit results for the iV° yield, mass and width with P-wave 
and S-wave assumptions. The fitted yields are found to be comparable with 
each other, while the mass and width show a systematic difference. We esti- 
mate systematic errors by varying the fitted yields by ±a for the background 
(±3) and B^ — * S°p7r+ (±11), and by taking into account the uncertainty in 
modeling the low mass (S°7r"'')x structure as discussed in the following. 

The simulated R{cos6p) distribution for B^ J]^N^ is almost fiat, as the 
generated cos 6p distribution is uniform for P- and S-waves. However, the 
R{cos9p) distribution is slightly affected by the assumed BW parameters due 
to efficiency changes in cos^p. We study the systematics of the fitted iV° mass 
and width due to the assumption on i?(cos^p), by changing the mass and the 
width in EvtGen in ranges between 1400MeV/c^ and 1620MeV/c^, and be- 
tween 200 MeV and 450 MeV, respectively. We find variations of ±lMeV/c^ 
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in the fitted mass and ±5 MeV/c^ in the width. We also study the systematic 
errors due to the parameterization of the (S^^tt"^) low mass structure. Instead 
of assuming a model with a single Breit-Wigner (S°7r"'')x, we consider a combi- 
nation of known states A* ^ Sc(2455)°7r+; A+(2625) (Ftotai < 1.9 MeV) [JS], 
A+(2765) ^ A+TT+TT- (r ~ 50MeV/c2), A+(2880) {T = 5.8 ± 1.3 MeV), and 
A+(2940) (r = IStfUeV). Here we use the partial widths for S°7r+ decay 
of the last three states given by Ref. [21]. We make a fit to the iV° mass, 
width and the yield of — > S°iV° with the individual A*+ yields floated and 
with the background fixed as mentioned previously. We obtain N{pTi^) mass 
and width values in good agreement with those obtained by the fit with the 
fiO ^ (S07r+)xp model. 

The branching fraction product B{B^ T,^N^)xB{N^ pir'^) is calculated 
as Ns/{Nb^ ■ e ■ CF ■ assuming Nb+b- = N^obo- For A^^ we use 

the P-wave yield in Table [2] as it gives a better confidence level than an S-wave 
fit. We use N^b = (387.7 ± 4.8) x 10^ for the integrated luminosity of 357 
fb"^, and the signal efficiency e = (5.18 ± 0.13)% from the MC simulation of 
5° S°iV°. We apply a correction factor CF = (86.7 ± 7.9)%, which takes 
into account the systematic difference in particle identification (PID) between 
data and MC simulation. Correction factors for proton, kaon and pion tracks 
are determined from a comparison of data and MC simulation for large samples 
of D*~^ D'^{Ktt)tt^ and A pTi~ decays. The overall PID correction factor 
is then calculated as a linear sum over the six tracks for the selected B signal 
events. We assign an error of 7.2% due to track reconstruction efficiency for 
the six charged tracks in the final state. The systematic error on the branching 
fraction arising from a quadratic sum of the uncertainties on N^bi the signal 
efficiency e, and particle identification CF and track reconstruction, is found 
to be 12%. Including the systematic error in the B^ —>■ Zl°iV° yield, we arrive 
at the total systematic uncertainty in the branching fraction of 17.6 %. Thus, 
we obtain the branching fraction product of B{B^ —>■ Sc(2455)°iV°) x B{N^ 
p7r+) = (0.80±0.15(stat.) ±0.14(syst.) ±0.21)xl0-^ and a significance of 6.1 
standard deviations including systematics. The last error is due to an uncer- 
tainty in i3(A+ ^ pK-Tr+) = (5.0 ± 1.3)% [H]. 

Next, we investigate goodness-of-fits with the masses and widths fixed to 
representative values for N{pTT^) states [18], and by floating the yields for 
B^ — > TP^N^ and B^ (S^7r+)xp. The fit results are summarized in Ta- 
ble [3l Here Iv2/,25 stands for a resonance of isospin / and spin S with an 
orbital angular momentum of L=S, P and D for L=0, 1 and 2, respectively. 
We exclude A states such as A(1600)P33 and A(1620)5'3i, as we have no sig- 
nificant structure in the pTT~ mass distribution in B^ S++p7r^ decay. The 
fits favor A^(1440)Pii and A^(1535)5n, while they disfavor A^(1520)L'i3 and 
A^(1650)5'ii. In the decay B^ — > S°iV° (assuming 5'(S°) = i), one expects 
a uniform cos 6'p distribution for the A^(1440)Pii state, A^(1535)S'ii state and 
A^(1650)5ii state, and a (1 + 3cos2^) distribution for the A^(1520)P'i3 state. 
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As shown in Fig. Mjo), the distribution has a peak only in the forward direc- 
tion, which is well reproduced by the mode (S°7r"'")xp. The remaining 
uniform distribution is due to S°iV°. Thus, the observed cos 6'p dis- 
tribution is consistent with both, iV(1440)Pii and iV(1535)S'ii states, with a 
preference for the former due to the width of the state. 

Table 3 

Results of the fits using the parameters of known N{pTT~^) resonances [18]. 



States 




Mass 


r 




(S°vr+)xP 








MeV/c2 


MeV 


Events 


Events 




N(1440) 


Pu 


1440 


300 


65± 10 


39 ±9 


37.6/34 


N(1520) 


Dvi 


1520 


115 


46 ±9 


53± 10 


53.5/34 


N(1535) 


5*11 


1535 


150 


58± 10 


43± 10 


40.1/34 


N(1650) 


5*11 


1655 


165 


44± 10 


55± 11 


74.2/34 



Finally, we try to perform a fit with an incoherent sum of the two Breit- 
Wigners, as we find that the fit results favor A^(1440)Pii and A^(1535)S'ii, 
and both give a distribution uniform in cos 6'p. Figure [7] shows the result of 
a fit to (a) the M{pTT^) and (b) cos 6'p distributions, where the N masses 
and widths are fixed to the values in Ref. |18j, and the individual yields are 
fioated. The histograms show the contributions from A^(1440) (solid), A^(1535) 
(dashed) states, B^ — > (S°7r+)xP (dotted), and the background (shaded). 
The yields are (37 ± 12) for the A^(1440) and (30 ± 11) for the iV(1535), 
while the B^ (S°7r"'")xp yield is (35 ± 9). We obtain the goodness of fit 
X^/ndf=30.3/33, which indicates a slight preference (by 2.7 a) for a mixed 
state of A^(1440) and Ar(1535) [18j. 

In summary, we study the three-body decay B^ Sc(2455)°p7r"'" with the 
same data set used for the analysis of the four-body decay B^ A+pvT+vr" [T]. 
We observe a broad pir^ mass structure near 1.5 GeV/ c^, and a uniform cos Op 
distribution with a sharp forward peak. To explain these structures, we find 
that contributions from an intermediate two-body decay B^ TP^N^ , non- 
resonant three-body decay B^ TP^pir^ and a low mass structure near thresh- 
old B^ {TP^TT^)y^p are needed. We perform a simultaneous fit to the M{p'K^) 
and cos 6p distributions with those three modes, and determine the yield and 
the relativistic Breit-Wigner parameters of the state for B^ —>■ Ti^N^. We 
obtain the branching fraction product of B{B^ Sc(2455)°A^°) x B{N^ 
p7r+) = (0.80 ± 0.15(stat.) ± 0.14(syst.) ± 0.21) x 10""^ with a signal significance 
of 6.1 standard deviations including systematics. The fitted mass and width 
are consistent with A^(1440)Pii and A''(1535)S'ii; both states also produce a 
uniform helicity distribution that is in good agreement with the data. The 
structure is also consistent with an interpretation in terms of an admixture of 
these two states. 
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Fig. 7. Simultaneous fit to the M(pTT~^) and cos 9p distributions with the A^(1440)Pii 
and A^(1535)S'ii Breit-Wigners. The histograms indicate the contributions from the 
Pii (solid) and 5ii (dashed) states, — >■ (E^7r"'")xP (dotted), and the background 
(shaded) . 
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